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1.  Objectives 


'^The  d^jectives  for  the  work  reported  here  are: 

V 

(1)  the  inclusion  and  evaluation  of  alternative  geometric  definitions  for 
channel  segment  inter-flowi 


(ii)  the  evaluation  of  the  legitimacy  of  routing  channel  segments 

separately  downstream,  to  better  model  in  and  out  of  bank  flows,' 


—  Ji* 


(lit)  the  establishment  of  a  working  Fulda  data  set  for  evaluation  of 
out-of-bank  conveyance  schemes ^ 

*  * 

(iv)  the  installation  of  SUN  Work  Station  and  transfer  of  HYMO  to  a  UNIX 
operating  system. 


These  objectives  represent  the  continuation  into  the  third  year  of  the 
project  stated  under  DAJA45-85-C-0011. 
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2.  Work  undertaken  in  the  6  months  period  to  December  1984 
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2,.  1  Out-of-bank  flood  conveyance 
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Following  the  March  1987  report,  work  has  continued  on  the  research  of 
improved  out-of-bank  flow  estimation  procedures.  As  suggested  in  the  last 
report,  work  on  the  redefinition  of  channel  geometry  to  include  turbulent 
exchange  between  segments  has  been  completed  -  section  2.2.  The  concept  of 
the  separation  of  segmented  flow  has  now  been  extended  to  the  flood  routing 
models,  to  better  model  the  discrete  paths  taken  by  out-of-bank  flow. 

2. 2  Redefinition  of  channel  geometry  to  incorporate  in-and-out 
of  bank  flow 


In  the  flow  rating  curves  generated  by  the  Manning  equation,  there  is  no 
provision  for  any  interaction  between  computations  undertaken  in  the 
specified  channel  sections,  i.e.  the  channel  cross  sectional  sub-areas 
(segments)  are  assumed  independent  in  terms  of  the  hydraulic  calculations. 
Significant  errors  can  occur  in  the  discharge  estimation,  which  is  of 
particular  importance  at  the  point  of  flood  plain  inundation. 

Chow  (1959)  suggested  that  redefinition  Of  the  area  and  wetted  perimeter 
terms  may  be  undertaken  to  provide  an  improved  discharge  estimate  based  on 
the  Manning  Formula.  The  four  methods  described  in  Table  2.1  have  been 
incorporated  into  the  CMPRC  subroutine  of  MILHY,  using  the  definition 
diagram,  figure 
2.1. 


The  trials  were  repeated  with  the  fictional  algorithm,  equation  2.1,  which 
assumes  friction  varies  as  a  function  of  stage,  removed: 

n'  »  n  -  0.0025  Araa/wected  perimeter  2.1 

The  results  are  shown  in  Tables  2.2  to  2.10.  Table  2.10  compares  results 
with  the  Manning  adjustment  (equation  2.1).  Note  that: 


flood  plane  width  :  channel  width  ratio 

ratio  of  flood  plain  :  channel  friction  'n'  values 


From  these  results  it  can  be  concluded  rhat:- 
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The  four  alternative  geometric  methods  do  make  2  significant  impact 
on  flow  rates-  Most  produce  lower  flow  rates  indicating  that  some 
allowance  for  interaction  in  the  flow  between  segments  has  bean 


The  flood  plain/channel  width  and  friction  ratios  have  an  impact  on 
flow  rates.  It  would  seem,  therefore,  that  the  optimum  method 
depends  on  the  individual  geometry  patterns - 

The  Manning  'n'  adjustment,  Table  2.10,  produces  flow  rates  varying 
up  to  40Z  from  computations  without  the  adjustment,  across  all 
geometries,  indicating  the  sensitivity  of  the  Manning  equation  to 


2.3  Application  to  hydrograoh  generation 


As  a  preliminary  study,  the  four  alternative  methods  have  each  been 
incorporated  into  HYM02  to  generate  a  hydrograph.  With  a  precipitation 
storm  of  2"  in  0.5  hours,  in  order  to  exceed  bankfull  conditions,  the 
hydrograph  was  routed  through  a  subcatchment,  using  each  of  the  fqur 
methods . 


Figure  2.2  shows  the  outflow  hvdrograph  for  HTM02  and  method  1,  whilst 
figure  2.3  Illustrates  methods  3  and  4,  where  in  both  cases  a  »  10,  b  »  4. 

When  the  flood  plain  channel  friction  ratios  were  at  their  greatest,  peak 
discharge  values  showed  a  variation  of  up  to  25X. 
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The  results  demonstrate  that  segment  turbulent  exchange  is  significant  and 
therefore  how  critical  the  rating  curve  generation  method  is,  especially  if 
flood  inundation  areas  are  to  be  predicted  with  acceptable  accuracy.  In 


3  “1 

figure  2.2,  a  discrepancy  of  10b  s  ,  indicates  an  elevation  difference  in  a 
typical  German  catchment  of  15  cm.  The  optimum  scheme  for  individual  valley 
geometry  and  accuracy  of  inundation  predictions  needs  to  be  investigated. 

For  the  present,  based  on  Knight  and  Hamed's  (1984)  results,  from  a  hardware 
application,  method  3  seems  the  most  appropriate  and  will  be  applied  here. 

2.4  Separate  routing  of  channel  and  flood  plain  water  . 

After  consultation  with  HALCROW  it  became  clear  there  is  a  discrepancy 
between  the  segmentation  of  water  in  the  rating  curve  generation  scheme  and 
the  lumping  of  all  water  in  the  routing  routine.  In  a  naturally  meandering 
channel,  there  is  a  tendency  for  flood  plain  water  to  "short-circuit"  along 
a  more  direct  route  (Fread,  1976).  Most  routing  models  treat  the  flood 
plain  as  stores  of  water  only,  or  as  in  the  HEC  model  generate  an  average 
routing  length  based  on  the  proportion  of  flow  on  the  flood  plains  and  in 
the  channel. 

At  the  present  time,  the  ROUTE  subroutine  is  being  investigated  in  order  to 
treat  each  segment  separately,  coincident  with  the  Chow  methods, 
acknowledging  the  turbulent  transfer  between  segments  illustrated  in  figure 
2.4.  The  PRTHY  subroutine  has  been  rewritten  enabling  it  to  recall  a 
previously  computed  rating  curve  and  thereby  convert  a  discharge  hydrograph 
to  a  stage  hydrograph. ,  A  series  of  trials  are  now  underway  to  pinpoint  an 

acceptable  discrepancy  in  stage  levels  between  segments  at  any  time 

interval.  It  is  envisaged  that  the  routing  method  will  involve  several 

"balancing  checks"  between  stations  to  prevent  the  build  up  of  water  in  any 

segment  at  any  time. 

Initial  results,  figure  2.5,  illustrate  the  necessity  of  this  investigation, 
where  the  outflow  hydrograph  for  the  same  reach  routed  with  the  segments 
together  (1)  or  divorced  (2),  are  dramatically  different.  Figure  2.6 
illustrates  the  same  reach,  routed  with  the  flood  plain  segments,  lengths 
reduced  by  5"  (1)  and  30"  (2).  Figures  2.7  and  2.8  show  a  different 
rainfall  storm,  routed  through  the  same  reach.  The  necessity  for  tfee 
separation  of  flood  plain  flows  are  therefore  clearly  dependent  on  the 
degree  of  "short  circuiting"  taking  place. 
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2.5  Establishment  of  The  Fulda  Data  Sec 

The  Fulda  River  data  set  has  been  established  after  a  visit  to  West  Germany 
earlier  this  year,  figure  2.9.  Data  was  collected  from  the  Water  Authority, 
Wasserwirtschaftsamt,  Fulda,  and  the  Meteorolgy  Office,  Deutcher 
Wetterdienst  Zentralant,  Offenbach,  Frankfurt. 

For  the  eight  gauging  stations  marked,  the  following  data  was  collected: 


-  cross-sectional  technical  drawing  of  gauging  site 

t 

-  stage/discharge  relationship 

-  stage  hydrographs  for  six  separate  storms,  recording  of  one  week  each 


During  the  visits,  the  gauging  sites  were  visited  and  valley  cross-sectional 
sketches  made,  valley  photographs  taken,  and  Hanning’s  ’n'  estimates  made. 

Meteorological  data  includes: 

-  daily  precipitation  values  for  the  period  three  weeks  proceeding  and 
the  week  of  the  stage  hydrograph  records 

-  continuously  recorded  precipitation  for  two  stations.  Bad  Hersfeld 
and  Kunzell-Dietershausen  for  five  of  the  six  storms  identified 
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daily  maximum  and  minimum  temperatures  and  snowfall  records  for  one 


Additional  data  alreadv  available  consists  of: 


topographic  maps 

soil  classification  maps 

flood  inundation  naps 


A  full  copy  of  this  raw  data  has  been  forwarded  to  Environmental  Laboratory, 

0 

W.E.S.,  Vicksburg.  Mississippi. 
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The  Fulda  data  set  is  now  running  with  the  completion  of:- 

1)  digitising  stage  hydrographs  to  15  minute  interval 

2)  parametric  regression  of  the  stage/discharge  relationship  and 
conversion  of  the  stage  hydrographs  to  discharge  values 

3)  combination  of  valley  and  gauge  site  cross-sections  with  topographic 
map  information  to  give  cross-sectional  data  for  direct  input  into 
HYMO 

4)  identification  of  sub-catchment  and  determination  of  the  channel 
network  in  a  logical  manner  for  the  routing  procedures 

5)  identification  of  the  proportional  major  soil  groups  for  each 
subcatchment 

6)  computation  of  weighted,  Thelsson,  precipitation  inputs  for  each 
subcatchment,  for  .each  of  the  identified  storms,  figure  2.10 

7)  construction  of  5  km,  radar  grid  and  preparation  of  precipitation  data 
for  moving  storm  analysis 

,  * 

3.  The  SUN  Work  Station 


A  considerable  advance  in  the  computing  application  of  HYMO  has  been  made 
with  the  transfer  of  all  work  to  a  SUN  station,  from  a  HONEYWELL  mainframe. 
The  3  MIPS  work  station  will  speed  the  development  of  individual 
subroutines,  and  the  numerous  applications  required  for  their  validation. 

Transfer  to  a  UNIX  operating  system  has  involved  minimal  changes  to  the 
inherited  Fortran  code,  and  HYMO  2  and  3  are  now  up  and  running.  The  full 
Fulda  data  set  has  also  been  transferred,  although  it  will  probably  be  two 
months  before  the  full  potential  of  this  machine  is  utilised. 
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4.  Research  plan  for  the  next  six  months 

(1)  Continuing  familiarization  with  the  SUN  work  station. 

(2)  Determination  of  the  acceptable  discrepancy  in  water  surface  elevation 
between  independently  routed  segments,  using  a  hypothetical  reach. 

This  will  be  used  to  establish  the  intermediate  reach  routing 
interval,  at  which  point  transfer  of  water  between  segments  will 


(3)  The  routing  (ROUTE)  and  travel  timetable  generation  (OiPTT) 
subroutines  will  be  rewritten  to  incorporate  the  option  of  discrete 
segment  routing. 

(4)  Sensitivity  analysis  of  ROUTE,  CMPTT  abd  CMPXC  subroutines,  combining 
the  segmented  routing  and  rating  curve  generation  schenes.  Thresholds 
in  stream  sinuosity,  flood  plain/channel  geometry,  and  frictional  ’n' 
values,  will  be  investigated  using  the  Fulda  catchment  data. 

It  is  anticipated  that  by  the  1  March  198S,  the  HYMO  model  will  incorporate: 

i)  the  five  alternative  cross-sectional  geometry  definitions  to 
incorporate  segment  turbulent  exchange 

ii)  discrete  segment  routing 

iii)  conversion  of  discharge  hydrographs  to  stage  hvdrographs 

The  most  appropriate  technique  (of,  for  example,  cross  sectional  geometry) 
will  wherever  possible  be  selected  by  the  program  using  the  thresholds  which 
will  have  been  previously  determined  in  the  sensitivity  analysis.  The 
operator  will  have  the  option  of  selection  but  will  be  expected  to  specify 
the  required  accuracy  of  output  data  required. 

C 

A  prototype  of  this  model  will  be  presented  during  a  visit  to  Vicksburg  in 
late  March  19SS. 
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Table  2.1 


Alternative  geometric  definitions  to  Incorporate  segment  intractions 

Chow  (1959) 
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Comparison  of  channel  geometry  redefir.iti 
methods  to  incorporate  inter-segment  flow 
h  -  HYM02  method 
i  -  Chow  method  1 
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Comparison  of  channel  geometry  redefinition 
methods  to  incorporate  inter-segment  flow 

3  -  Chow  method  3 

4  -  Chow  method  4 
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